Numerical Analyses on 
Moduli Space of Vacua 

<N 

^ Jonathan Hauenstein^ Yang-Hui He^ & Dhagash Mehta^ 

<N 

o 

^ Department of Mathematics, North Carolina State University, 
^ Raleigh, NC 27695-8205, USA 
hauenstein@ncsu.edu 

+Y» ^ School of Physics, NanKai University, Tianjin, 300071, P.R. China; 

Department of Mathematics, City University, London, ECl V OHB, UK; 
^ (— I Merton College, University of Oxford, OXI4JD, UK 

hey ©maths, ox. ac. uk 

K*" ^ Department of Physics, Syracuse University, 

00 

cn Syracuse, NY 13244, USA. 



O 



(N 
^ — I 

> 
>< 



dbmehta@syr. edu 



Abstract 

We propose a new computational method to understand the vacuum moduh space of 
(supersymmetric) field theories. By combining numerical algebraic geometry (NAG) 
and elimination theory, we develop a powerful, efficient, and parallelizable algorithm to 
extract important information such as the dimension, branch structure, Hilbert series 
and subsequent operator counting, as well as variation according to coupling constants 
and mass parameters. We illustrate this method on a host of examples from gauge 
theory, string theory, and algebraic geometry. 
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1 Introduction 



The vacuum is undoubtedly one of the key ingredients in any quantum field theory, 
therefrom arises particles and dynamics. Vacua which have non-trivial geometry are 
familiar to us, e.g., the famous "Mexican hat" of the Higgs potential has its true 
minimum at a circle rather than a mere point (a local minimum at the origin). In 
general, the expectation values of scalar fields can parameterize intricate manifolds 
called the vacuum moduli space. In field theories with supersymmetry, where such 
scalar fields abound, the situation become particularly pronounced and complicated 
geometries can arise as the vacuum moduli space. In the context of string theory, 
realizations and interpretations of the geometry of the vacuum as the low-energy limit 
of a compactification or holographic scenario is key to the engineering of field theories. 
Indeed, one could even begin to unravel unexpected structure in the standard model 
by examining the vacuum under this light [l}]3]. Thus, to have a mathematical tool to 
study the geometry, ranging from topological to metrical issues, of the moduli space of 
vacua of field theories is therefore much needed. 

Computational algebraic geometry has become one of the most useful tools to study 
numerous phenomena in theoretical physics. In the recent years, active research has 
been concentrated on the rich interplay between algorithmic geometry and theoretical 
physics, especially in gauge and string theory [4j. With the increasingly powerful 
computers the methods and algorithms to use algebraic geometry in practice to solve 
problems arising from theoretical physics have been becoming more relevant than ever. 
Harnessing these developments to study the moduli space of vacua has been initiated 
over the past few years l2l[5]. 



The traditional computational methods are based on symbolic computational alge- 
braic geometry, most of whose sub-methods and sub-algorithms rely on Grobner basis 
techniques (cf. j|6] for a nice introduction for physicists). Roughly speaking, given a 
vanishing set of the polynomials, the so-called Buchberger Algorithm (BA) or its refined 
variants compute a new equivalent system of polynomials, called a Grobner basis [T], 
which has nicer properties; this is analogous to Gaussian elimination for linear systems. 
Nowadays, efficient variants of the BA are available, e.g., F4 |8], F5 p|, and Involution 



Algorithms 10 . Symbolic computation packages such as Mathematica, Maple, Reduce, 



etc., have built-in commands to calculate a Grobner basis. Moreover, Singular 11 



COCOA 12 , and Macaulay2 13 are specialized packages for computational algebraic 



geometry, available as freeware, and MAGMA 14 is also such a specialized package 
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available commercially. 



In 15 -17 , the Grobner basis method was used to answering questions pertinent 
to string and particle phenomenology and a publicly available computational package, 
called Stringvacua, was designed to interface Mathematica with Singular for convenient 
usage. Utlizing Stringvacua, extracting such important information as the dimension 
of the vacuum, the number of real roots in the system, stability and supersymmetry 
of the potential, or the branches of moduli space of vacua, etc. can be readily settled 
using only a regular desktop machine in many circumstances. 

Nonetheless, there are a few well-known problems with such methods. First, the 
BA is known to suffer from exponential complexity: the computation time and the 
RAM required by the machine increases exponentially with the number of variables, 
equations, degree, and terms in each polynomial. Moreover, BA is usually less efficient 
for systems with irrational coefficients. One habitually has to resort to randomizing 
over the space of integer (or prime) coefficients and work over finite fields. Another 
shortcoming is the highly sequential nature of the BA. 

Although Grobner methods have been an immensely useful tool to study various 
theoretical physics problems, more prominently so in string and gauge theories, the 
time has come when one needs to depart from the computation of relatively simple 
models and develop approaches for attacking more elaborate systems. Unfortunately, 
methods based on BA often run out of the steam for more realistic models due to the 
aforementioned shortcomings. 

A recently developed approach called numerical algebraic geometry (NAG) over- 
comes many of the shortcomings of Grobner basis techniques. The core method in NAG 
is numerical polynomial homotopy continuation (NPHC): given a system of polynomial 
equations which is known, one computes its isolated solutions by first solving a related 
system and then tracking "solution paths". More specifically, one first estimates an 
upper bound on the number of isolated solutions of the given system, creates a new 
system which is easy to solve and has the same number of isolated solutions as this 
estimated upper bound, and finally one tracks the paths from each of the solutions of 
the new system to the original system. The paths which converge yield solutions of the 
original system. In this way, for a system which is known to have only finitely many 
solutions, one can obtain all its solutions. This is a remarkable method in its own right 
based on this guarantee to find all of the solutions unlike other numerical methods such 
as the Newton-Raphson or its sophisticated variants. Unlike standard Grobner basis 
techniques, the NPHC method is embarrassingly parallelizeable and hence one can often 
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solve more complicated systems efficiently using computer clusters. 

In many theoretical physics problems, one encounters systems of equations that 
have infinitely many solutions. In this case, the problem of solving such a system is 
transformed to a collection of systems of equations which have isolated solutions corre- 
sponding to the components of the original system and then invoke the NPHC method. 
There are several sophisticated numerical packages well-equipped with path trackers 
such as Bertini [l8|, PHCpack (l9), PHoM f20], HOMPACK |2l), and HOM4PS2 |22||23 
which are all available as freewares from the respective research groups. 



NAG was introduced in particle theory and statistical mechanics areas in Ref. 24 



Subsequently, the NPHC method was used to solve systems arising in numerous phys- 
ical phenomena in lattice field theories 25 26 , statistical physics 27 30 , particle 



phenomenology 31 , and string phenomenology 32 



Recently, we undertook the systematic introduction of efficient methods in NAG 
to the investigation of classes of problems in gauge and string theories p53]. There we 
focused particularly on systems of equations with finitely many solutions as these are 
important to phenomenological questions such as looking for isolated extrema of poten- 
tials. In some sense, what we undertake in our present study of moduli space of vacua 
is the substantial generalization thereof where we can proceed from zero-dimensional 
solution sets to solution sets of arbitrary dimension. We show that studying moduli 
spaces from an algebraic geometric perspective reduces to a problem in elimination 
theory and then utilize the power of NAG to address this problem. 

In particular, we propose a complete algorithm in this paper for computing the 
irreducible decomposition of an elimination ideal, which is based on recently developed 
techniques regarding projections in NAG 34,35 . We also combine classical methods 
in algebraic geometry with NAG to compute the Hilbert function and Hilbert series of 
components of the solution set which satisfy a certain regularity condition that holds 
for a wide range of problems in quantum field theory, namely the components are 
arithmetically Cohen-Macaulay (see p6 , Chap. 1] for more details). This is a mild 
technical constraint which most algebraic varieties we encounter in physics obey. 

The outline of the paper is as follows. We begin, in Section |2} with a pedagogical 
introduction of vacuum moduli spaces in supersymmetric field theories and explain why 
elimination theory of ideals is the most efficient approach. Next, we turn to the mathe- 
matical problem of numerical elimination theory and propose a novel algorithm to com- 
pute the irreducible decomposition of elimination ideals in Section |3j In Section |4} we 
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return to the physics and attack various exphcit problems using our algorithm, ranging 
from quiver gauge theories to Calabi-Yau manifolds and from string phenomenology to 
instantons. These demonstrate that our numerical elimination algorithm is a powerful 
approach by solving problems well beyond the reach of traditional methods. 



2 Supersymmetric Vacuum Moduli Space 

As emphasized in the introduction, the study of the moduli space of vacua for gauge the- 
ories, especially ones of supersymmetry in (3+l)-dimensions, is of crucial importance 
to such diverse disciplines as string theory, particle phenomenology and mathematics. 
In this section, we will therefore turn to a fairly pedagogical presentation of computing 
the vacuum moduli space, reflecting the three paradigmatic shifts in approach: the 



original problem (cf. [37J), the symplectic quotient method of Luty- Taylor 38 , and 
the algebro-geometric algorithm of [Tlls]. 



Our starting point is the action for an A/" = 1 globally supersymmetric gauge 
theory in 3 + 1-dimensions, written as an integral in superspace (with repeated indices 
implicitly summed): 



s = j d^x j d^e $Je^$, + (^-^ j d^e ttWo^w + j d^e + h.c. 



(2.1) 



Here, the $j are chiral superfields transforming in some representation Ri of gauge 
group G] \^ is a vector superfield transforming in the Lie algebra g of G; Wa = 
iD e~^Dae^ , the gauge field strength, is a chiral spinor superfield; and iy($) is the 
superpotential, which is a holomorphic function of the 

The supersymmetric vacuum moduli space (VMS) Ai is simply the set of field 
configurations which minimize the effective potential obtained from integrating out 



the Grassman variables 9i in (|2.1|). Explicitly [39|, they are the solutions to 

JL 



F-Flatness : = {^W{(pi) = 0}; 



D-Flatnes : = = • ^^'^^ 

i 

with (pi being the vacuum expectation values of the scalar components of $j. In the D- 
flatness conditions, we have chosen the Wess-Zumino gauge and are the generators 
of Q. The index i runs over the fields and A, over the gauge group generators. 



There is gauge redundancy to the action in (2.1), namely $ — )■ exp(iA) ■ $; e — )■ 



exp(— zA^) ■ ■ exp(— iA), in terms of the parameter A G 0. Choosing A to be complex 



makes the full gauge symmetry to be the complexification of the original gauge group 
G. Following a series of systematic studies of fixing this gauge symmetry 40 43 , it 



was realized in 38 that for every solution to the F-flatness constraints there is one 



and only one solution to the D-flatness conditions. Therefore, D-flatness is merely 



a gauge-fixing condition and the full set of solutions to (2.2) is simply a symplectic 
quotient of the space of F-flatness by the complexified group G"^. In other words, as an 
algebraic variety, the VMS is the GIT quotient 

M=T II G\ (2.3) 

From a computational algebraic geometric point of view, this rather abstract quo- 
tient was made explicit in j^. Let there be n superfields $j=i,...,n with associated vac- 
uum expectation value of scalar component 0j and polynomial superpotential W{(^i). 
This allows us to define the F-fiatness equations as the Jacobian ideal (^W^) in the 
polynomial ring i? = C[</)i, . . . , 

Since the D-fiatness equations are gauge fixing conditions, we can consider the 
set of gauge invariants, composed of polynomials in 0j. These are traces of matrix 
products so that the final answer will carry no free gauge index. Moreover, there 
should be a minimal generating set of the ring of such invariants, which we denote as 
D = {rj{{(j)i})} where rj=i^...^fc are polynomials and k is a positive integer (potential 
quite large). Hence, we have another polynomial ring S = C[ri, . . . , r^]. 

The realization in \2j is that the polynomial map D is a ring map from the quotient 
ring ~ i? y (^^) to ™S ^■ 

lid \ D={r,{{^,})} 

C[0i,...,0„] /(— ly) y C[ri,...,rfc] . (2.4) 



The image of this ring map is then the VMS, defined as an affine variety in S: 

M ~ Im(J- A S) . (2.5) 



The above procedure is perfectly adapted for the language of computational com- 
mutative algebra, especially using either of the wonderful computer packages Macau lay 



2 [13] and Singular 44 . However, the Grobner basis algorithm is central to this com- 
putation and is precisely the computationally intensive step we wish to circumvent. To 
this end, one can use a standard trick in polynomial manipulations to re-phrase and 
simplify the above process. Instead of considering D as a map from R to 5, we can 
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simply enlarge the ring R by adjoining new variables ?/j=i,...,fc to form the larger ring R 
and consider the following ideal therein: 

/ dW \ 

-^i > C = C[0i=i„„,„,?/j=i,...,fc] . (2.6) 



The image in (2.5) describes the relations (first syzygy) amongst the variables Vj 
subject to the F-flatness conditions of J-". This is equivalent to systematic elimination 
of the 0j variables in R so that only the Uj variables remain. The resulting ideal in R 
will explicitly have only polynomial relations amongst the Uj variables. Since we have 
set each yj equal to ^"^({04}), these relations define precisely the VMS in R. Hence, this 
formulation bypasses the language of ring maps entirely and we are now confronted 
with an elimination problem. 

Indeed, there are classical and modern techniques for elimination, such as using 
resultants and computing the Grobner basis of the ideal with a so-called elimination 
ordering for the variables. Although it is impossible for us to provide a complete 
account of techniques in elimination theory, we will provide a little historical perspec- 
tive. Gaussian elimination is the standard elimination technique for systems of linear 
equations. In simple terms, Buchberger's algorithm can be thought of as a generaliza- 
tion of Gaussian elimination for polynomial systems. As mentioned above, resultants 
were classically applied to perform eliminations and there are many renowned math- 
ematicians associated with them including Bezout, Euler, Sylvester, Cayley, Dixon, 
and Macaulay, and reinvigorated by hazard. Another recent approach that is related 
to, but independent of, Grobner bases is Wu's characteristic sets. We refer the reader 



to 45 and the references therein for more details on these methods and related history. 



In this paper, continuing with the experience gained from 33 , we will try to avoid 
the expensive step of computing Grobner bases when all the information needed for the 
VMS are such quantities as the dimension, irreducible components, and even Hilbert 
series. Instead, we will harness the power of numerical algebraic geometry in § |3] to 
perform this computation. 

2.1 Summary of Algorithm 

It is expedient to summarize the algorithm which we employ to compute the VMS: 
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• INPUT: 

1. Superpotential W{{(f)i}), a polynomial in variables 0i=i,,,.,„, 

2. Generators of gauge invaraints: rj{(f)i), j = 1, . . . , k polynomials in (pf, 

• ALGORITHM: 

1. Define the polynomial ring R = C[(f)i=i^,„^n,yj=i,...,k], 

2. Consider the ideal I = (|^; yj — rj{(j)i)), 

3. Eliminate all variables 0j from I G R, giving the ideal M in terms of yj; 

• OUTPUT: M corresponds to the VMS as an affine variety in C[j/i, . . . , 7/^]. 



3 Numerical Algebraic Geometry and Elimination 

Having phrased our physics as a problem in mathematics, in this section, we will 
develop the necessary mathematical and computational tools. The Grobner basis ap- 
proach presented in § [2] uses techniques from algebra via manipulation of equations to 
compute eliminations. The numerical algebraic geometric techniques presented here 
compute eliminations based on geometry via points. The first two techniques presented 
are for computing low degree polynomials, which can be used to improve both Grobner 
basis and additional numerical algebraic geometric computations. The other two tech- 



niques utilize the numerical algebraic geometric notion of witness sets (see 46 , Chap. 



13]) to compute the elimination. We close the section by computing the Hilbert series 



in the arithmetically Cohen-Macaulay case (for more details, see 36, Chap. 1]). An 



ongoing project is a practical method for computing the Hilbert series in general. 

For a polynomial system F, define V{F) = {x \ F{x) = 0} to be the variety 
defined by F. To simplify the presentation, throughout this section, we let 7r(x, y) = y 
and Z = V{Q) for some polynomial system Q. We aim to compute the elimination 
Y = 71 (Z) which will be considered based on the following three cases: 

1. Parameterization: Q{x,y) = y — F{x); 



2. Restricted parameterization: Q{x,y) = 

3. General case: Q{x,y) is some polynomial system 



G{x) 

y - F{x) 



and 
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where F and G are polynomial systems. In the parameterized case, Z is an irreducible 
and smooth variety implying that Y is also irreducible. In the other two cases, irre- 
ducibility depends upon the structure of the polynomial system Q. Additionally, Q 
may impose multiplicities on the irreducible components of Z, e.g., Q{x,y) = (x — y)^. 
However, by using isosingular deflation 47 where appropriate, we can assume without 
loss of generality that each irreducible component of Z has multiplicity one with re- 
spect to Q. This allows us to compute dimensions simply by computing the dimension 
of a linear space related to the Jacobian of Q as shown in the following section. 



3.1 Dimension of the image 



Let V G Z = V{Q) he an irreducible component of multiplicity one with respect to Q 
and (x, y) G x C^"'^ be a general point of V. Then, the dimension of the irreducible 
algebraic set vr(V^) (which may or may not be an irreducible component of Y = vr(Z)) 
is easily computed via linear algebra ^| as follows: 

dim7r(V^) = dimnullJ^(x, y) — dimnullJ^(x, _i;7v] (3.1) 

where JQ{x,y) is the Jacobian matrix of Q at {x,y) and JQ{x,y)[... ^i:]\f] is the matrix 
consisting of the first columns of JQ{x,y). In the parameterized case, namely 
Q{x,y) = y — F{x), we know that Z and Y are both irreducible with dimZ = and 

dimF = A^ - dimnullJF(i). (3.2) 



EXAMPLE 1 We consider the parameterized example Q{x,y) = y — F{x) where 



F{xi, X2, X3, X4, X5, Xq, X7, Xg) 



X1X5 + 

XsXs + X^Xq 
XiXj + 0:2X8 

X3X7 + X4X8 
X6X7 - X5X8 



The set Z = V{Q) is irreducible of dimension 8. The dimension of Y 
8 — 3 = 5 since dimnullJ-F(x) = 3 for a random x G C^. 



7r(Z) is 
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3.2 Low degree polynomials via interpolation 



In order to use interpolation techniques to compute low degree polynomials, one must 
have the ability to compute points on Y. In the parameterized case, computing points 
on Y is trivial since F{x) G Y for every x. Under the assumption that Z is irreducible, 



the numerical algebraic geometric technique of sampling an irreducible component 46 
§15.2] allows one to compute arbitrarily many points on Y given a sufficiently general 
point {x, y) & Z. li Z is not irreducible, then we need a sufficiently general point on 
each irreducible component of Z. In short, sampling Z, and hence Y via the projection 
TT, reduces to following a curve starting at (x, y) defined by the intersection of Z with 
a family of linear spaces of complimentary dimension to Z. Additionally, we restrict 
our attention to low degree polynomials due to the exponential growth of the matrices 
involved and to avoid issues regarding numerical stability. Nonetheless, computing 
some low degree polynomials can be used to expedite further computations. 

Given a finite dimensional linear space of polynomials V, interpolation using suffi- 
ciently many points on each of the irreducible components of Y yields the subspace of 
polynomials in V which vanish on Y. This is accomplished by computing the null space 
of a matrix constructed by evaluating the points at a basis for V. We note that using 
more points, thus adding additional rows to this matrix, and rescaling each row will 
often drastically improve the conditioning of this matrix, e.g., see ^8]. Additionally, 
once such a polynomial is computed, it is often easy to use either symbolic evaluation 
or Grobner basis methods to validate that a given polynomial does indeed vanish on Y. 

EXAMPLE 2 For the elimination problem in Ex. [I| the set Y was verified to be a 
hypersurface in meaning that Y = V(/) for some polynomial f . Before we consider 
using interpolation to compute this polynomial f , we note that since each polynomial 
in F is homogeneous of degree 2, f is also homogeneous. Thus, we will restrict our 
attention to homogeneous polynomials in yi, . . . ,yQ. 

We first consider V to be the set of linear homogeneous polynomials in yi, . . . ,yQ 
taking the variables as the basis. Each row of the matrix under consideration is thus 
of the form 

F,{x) F^ix) F,{x) F,{x) F,{x) F,{x) ] (3.3) 

where x G since yi = Fi{x). After picking 6 random points in C*, it is easy to 
verify that the resulting 6x6 matrix has full rank meaning that there are no linear 
polynomials which vanish on Y . 
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We then consider V to be the set of homogeneous polynomials in yi, . . . ,yQ of degree 



2 taking the monomials of degree 2, namely yl,yiy2,yiy3, ^ 
row of the matrix under consider is thus of the form 



as the basis. Each 



Fg(x) 



(3.4) 



where a; G C^. After picking 21 random point in C^, it is easy to verify that the resulting 
21 X 21 matrix has rank 20 meaning that there is a quadratic polynomial which vanishes 
on Y . Computing the null space reveals this polynomial to be yiy^ — 7/2^5 + ysVi- By 
substitution, it is easy to verify that FiFq — F2F5 + F3F4 = 0. 



3.3 Low degree polynomials via lattice base reduction 

If the polynomials defining Z have rational coefficients, then the same is true for Y. 
The immediately follows from the fact that methods involving only Grobner bases 
maintain the same field of coefficients. As described in (49| § 3.3], lattice base reduc- 



tion techniques, such as LLL 50 and PSLQ 51 , provide an approach to compute 
polynomials vanishing on Y. For this technique, only one sufficiently general point is 
needed for each irreducible component of Y. 



Fix a sufficiently general point y on the irreducible component V CY. As in § |3.2 
we want to compute the subspace of polynomials in a finite dimensional linear space of 
polynomials V which vanish on V. Let v be the vector obtained by evaluating y on a. 
basis for V. Lattice base reduction techniques compute an integer vector w such that 
f ■ w ~ 0. By changing the precision and recomputing, one is able to reliably determine 
if w actually corresponds to a polynomial which vanishes on V. Moreover, as stated 
above, validation is often trivial. 

EXAMPLE 3 Consider repeating the computation from Ex.^using lattice base reduc- 
tion rather than numerical interpolation. We proceed by first fixing a random x G C^. 



For the linear polynomials, we apply PSLQ to the vector computed in (3.3). Us- 
ing 20 digits, this yields w = [—431,-156,571,300,597,95] which changes to w = 
[8994,16067,-12869,-33044,27787,56416] when using 30 digits. This suggests that 
there are no linear relations with small integer coefficients. 



For the quadratic polynomials, PSLQ applied to the vector computed in (3.4), using 



both 20 and 30 digits, computes the vector w corresponding to yiy^ — 2/21/5 + ysy^. 
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3.4 Solve by slicing 



Suppose that we are in the restricted parameterization case, which includes the pa- 
rameterization case by simply taking G{x) to be the zero polynomial. We assume 
G : — 7- C" and F : — )■ C*^ are polynomial systems and V C V{G) is an irre- 
ducible component (where V = in the parameterized case). By using isosingular 



deflation 47 if needed, we can assume that V has multiplicity one with respect to G. 



The key numerical data structure for y is a witness set which is the triple {G, L, W} 
where L : — )■ is a system of general linear polynomials where i = dimV^ and 
W = V n V{L). The set W, called a witness point set, consists of deg V" points. 



The set F(V) is irreducible since V is irreducible and F is polynomial. Its dimension 



can be computed using (3.1), say i' = dimF(l^). Since we do not have a system of 



polynomials for which F{V) is an irreducible component, we can not readily compute 
a witness set for F{V). However, we are able to compute a pseudo-witness set l34j for 



F(y) which allows us to answer questions about F{V), such as its degree. We do this 
by considering the following polynomial system constructed by adding slices to Q: 



n{x,y) 



G{x) 
C{x) 

y - F{x) 



(3.5) 



where £ : ^ C^-^' and C : ^ C^' are systems of general linear polynomials. 
Let W be the set consisting of the finitely many solutions of "H whose x coordinates 
lie on V. Such a set may be computed efficiently from a witness set for V using 

C(x) 



regeneration 



52 



. The quadruple {Q, it, M, W} where M(x, y) 



is a pseudo- 



witness set for F{V) with degF{V) = |7r(W^)|. For any point y G F{V), the fiber over 
?/ in is the set 

Fy = {xeV\y = F{x)}. 
For a general point y e F{V), dimFj^ = £ - £' and degFj^ = |M/|/|7r(iy)|. 

EXAMPLE 4 For the elimination problem in Ex. Q the set Y = F(C8) = 7i{Z) was 
verified to be a hypersurface in so that i = 8, i' = 5, and i — i' = 3. To compute 
the degree of Y and the degree of the general fiber, we solve 'H{x,y) = defined in 



(3.5) where we take G{x) to the zero polynomial. That is, Ti is a polynomial system 
depending on 14 variables that consists of 3 linear slices C, 6 polynomials y — F{x), 
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and 5 linear slices C The set W consists of A points with \7r{W)\ = 2. Therefore, 
degY = 2 and the general fiber has dimension 3 and degree 4/2 = 2. 

To compute the irreducible components of Y, we first compute a witness set for each 
irreducible component of V{G), called a numerical irreducible decomposition (see |46| 
Chap. 15] for more details). Then, for each irreducible component V C V{G), we 
follow the aforementioned approach to compute a pseudo-witness set for F{V). The 
irreducible components of Y are precisely the inclusion maximal sets of 

V C V(G) is an irreducible component | . 

These inclusion maximal sets can be computed using the pseudo-witness sets con- 



structed above with the membership test provided in 35 



3.5 General approach 



The general case follows a similar approach as above except that we must start with 
a numerical irreducible decomposition of Z = V{Q), that is, compute a witness set for 
each irreducible component of Z. For each irreducible component V, we follow the 
basic approach 



34 



to construct a pseudo-witness set for t^{V). Then, the irreducible 



components of Y = 7r(Z) are precisely the inclusion maximal sets of 



niV) 



V C V{G) is an irreducible component | 



which can be computed using the membership test of 35 



EXAMPLE 5 For a basic example, consider computing the closure of the set of all 
(a, c) G such that ax^ + 2ax + c has at least one root where the derivative with respect 
to X vanishes. That is, we want to compute Y 

ax"^ + 2ax + c 



g{x,a,c) 



n(Z) where 
V{Q), and 7!'{x,a,c) 



(a,c). 



2ax + 2a 

A numerical irreducible decomposition shows that Z decomposes into a line and quadratic 
curve, namely V\ = {(x, 0,0) \ x E C} and V2 = {(x,a, 1) | ax = —I}, respectively. 
Using 



3.1 



we have that dim7r(\4) = and dim7r(V2) = 1. Since 7r(Vi) is the sin- 
gleton containing the origin, a witness set for Vi trivially yields a pseudo-witness set 
for 7r{Vi). A pseudo-witness set for n {¥2) using [341, which simply moves linear slices, 
yields that deg7r(V2) = 1 and the general fiber consists of one point. In particular, 
7r(V2) = {(a, 1) I a G C}. The final step uses [351 to determine that tt{Vi) is not con- 
tained in 7r(V2) meaning Y has two irreducible components, namely a point and a line. 
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3.6 Hilbert series 



Positive dimensional arithmetically Cohen-Macaulay (aCM) schemes retain informa- 
tion under slicing by hyperplanes and hypersurfaces [36| Chap. 1]. In short, this implies 
that the Hilbert function and Hilbert series for an aCM scheme can be computed di- 
rectly from a witness point set. 

We divert briefly into projective algebraic geometry. Assume that V C P" is a 
closed subscheme of dimension r > 1. As above, we will assume that V is generically 
reduced, i.e., simply consider it as an algebraic set in P". There are two equivalent 
ways that we will describe that one can use to determine if V is aCM. The first, which 
is typically taken to be the definition of aCM, is that the quotient S/Iy is a Cohen- 
Macaulay ring where S is the homogeneous polynomial ring C[xo, . . . , x„] and /y is 
the saturated ideal associated to V. That is, the KruU dimension of S/ ly is the same 
as the depth of S/Iy- The second is that deficiency modules of V are trivial, namely 
{M'''){V) = Hi{Iv) = for 1 < z < r where Jy is the ideal sheaf of V. 

We denote the Hilbert function of V using an integer vector notation, namely 
Hy = {hy, hy, . . . ), and the Hilbert series as the function HSvif) = YlJLo ^y^"'- 

Let W he a witness point set for V, that is, there exists a general codimension r 
linear space C such that W = VnC. The Hilbert function for W, considered as simply 
a zero-dimensional set, can be computed via numerical linear algebra using Veronese 
embeddings of the points of W, e.g., see [is]. We know h^r = 1, h^y > for all j > 0, 
and Yl'jLo^w ~ 1^1 iinplying h^y = for all j sufficiently large. In particular, the 
Hilbert series for W, namely HSw{t), is a polynomial. 



The following theorem, which follows from Corollary 1.3.8(c) and Section 1.4 of 36 
shows how to compute the Hilbert series and Hilbert function for V from Hyy. 

THEOREM 3.5 Let V be a generically reduced and aCM scheme of dimension r > 0. 
Let W he a witness point set for V with Hilbert function Hyy = {h^r, h\y, . . . ). Then, 

k jl jr-l 

1. Hv = (/ly, /iy, . . . ) where hy = ^^- ■ - ^ h^w>' ^'^^ 

jl=0j2=0 jr=0 

2. HSvit) = 

Even though we have stated this for witness point sets of projective varieties, we can 
simply restrict to an affine patch by taking linear slices that are general with respect 
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to the chosen patch. In short, the same computation holds for affine varieties and a 
generahzation to pseudo-witness point sets is trivial. 



EXAMPLE 6 Reconsider the elimination problem first considered in Ex. [I| We know 
that Y is a hypersurface in of degree 2 that one can verify is aCM. Let £i, . . . , £5 
be general hyperplanes in and consider the witness point set = fl £1 fl ■ ■ ■ fl £5 
consisting of two points, say W = {wi,W2}- Since the matrix [wi W2] has rank 2, we 
know h^r + hly = 2 yielding Hw = (1, 1, 0, 0, ... ). Table^presents the Hilbert function 
and Hilbert series for Vk = V nCiH- ■ - nCk for k = 0, . . . ,5 which can be easily verified 
using either Macaulay 2 75' or Singular 44^- ^''^ particular, HSv{t) = (1 + t)/(l — t)^. 



k 



HS 



5 (1,1,0,0,0, 



1 + t 



:i,2,2,2,2. 



1+t 
i-t 



3 (1,3,5,7,9, 



1+t 



2 (1,4,9,16,25,... 



1+t 



1 (1,5,14,30,55,... 



1+t 



(1,6,20,50,105, 



1+t 



Table 1: Summary of Hilbert functions and Hilbert series 



4 Illustrative Examples and Applications 

Having outlined the method of attack in the computation of the moduli space of vacua 
and discussed in detail the mathematical background underlying the algorithm, we 
now, in this section, turn to concrete examples which arise in actual physical situa- 
tions and demonstrate how our numerical perspective affords an efficient outlook on 
understanding the physics. 



4.1 Warmup: an Orbifold Theory 

Let us start with a four-dimensional N" = 1 theory which is famous to the AdS / CFT 
literature. This is the world- volume theory of a D3-brane transverse to the non-compact 
Calabi-Yau threefold locally realized as the quotient of C'^ by Z3 which acts on the three 
complex coordinates {x,y,z) of the former by {x,y,z) — )■ exp{^){x,y, z) (otherwise 
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known as (1,1,1) action). The theory is given by a U{1)^ quiver together with a 
superpotential W: 




w= E 

a,/3,7=l 



via) yyfj) y- 



^23 



(7) 
31 ■ 



(4.1) 



There are nine fields xif^X',^^', and X. 



'-23 5 



1, 2, 3, with the subscript ij signi- 
fying an arrow from node i to j and superscript a denoting that there is a multiphcity 
of three arrows for each pair of nodes. Subsequentfy, there are 3^ = 27 gauge invariant 
operators formed by their products corresponding to the closed cycles in the quiver. In 
the superpotential, eap-y is the standard totally antisymmetric Levi-Civita symbol on 



three indices. Taking the partial derivatives of W with respect to the fields, we obtain 

3 3 3 

^(7) ^ vi'^) vi'y) o-K,^ ^ vio') yW 

"31 ' to/37^12 
a,7=l 



ea/37^23 



K7) 
"31 ' 



and ^ 6^/37X^2 ''-^23 



9 F-terms: 

/3,7=1 0:7=1 o,/3=l 

Of course the moduli space of vacua of F-fiatness and D-fiatness should give the 
affine equations for C^/Zs as an algebraic variety; this was studied in detail in §2.3.2 
of [2]. Here, in light of our new methodology, we are trying to solve the following 
problem (below, a, /3, 7 = 1, 2, 3 with free indices ranging fully): 



Given ideal 



Eliminate 



3 



^12 ^23 - 



'(7) 
"31 5 



fco/37-^23 ^31 ' 



3,7=1 
3 



-(«) x-{7) 



fco/37-^12 ^31 ' 



a, 7=1 
3 
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y{*) -V-{*) \^(*) 
^12 5 ^23 1 ^31 



(*) \-(*) 



(4.2) 



This ideal consists of 3^ + 3 ■ 3 = 36 generators in the polynomial ring with 3^ + 9 = 36 
variables. After elimination, we obtain relations amongst the 27 variables Da^-^. 

Our first step is to compute the numerical irreducible decomposition of the 9 gen- 
erators involving only the X variables, which yields one irreducible component of di- 



mension 5 and degree 6. Therefore, the eliminant is irreducible and, using (3.1), has 



dimension 3 with the general fiber having dimension 2. Using the approach of § 3.2 
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one readily obtains the following 17 linear and 27 quadratic equations: 



yii2 - 2/211 


= 2/113 - 2/311 


= 2/121 


- 2/211 = 


2/122 — 2/221 


— 2/123 ~ 2/321 






= 2/131 - 2/311 


= 2/132 


- 2/321 = 


2/133 - 2/331 


= 2/212 - 2/221 








— 2/223 


— 2/322 — 










= 2/233 — 2/332 


= 2/312 


- 2/321 = 


2/313 — 2/331 


= 2/323 - 2/332 = 0, 


2/332 


— 2/3222/333 


= 2/3312/332 — 


2/3212/333 


= 2/3222/332 


— 2/2222/333 


= 2/3212/332 — y22\yi,-iZ = 


2/331 


- 2/3112/333 


= 2/3112/332 - 


2/2112/333 


= 2/3222/331 


— 2/2212/333 


= 2/3212/331 - 2/2112/333 = 


2/322 


— 2/2222/332 


= 2/3112/331 - 


2/1112/333 


= 2/2222/331 


— 2/2212/332 


= 2/2212/331 — 2/2112/332 = 


2/321 


~ 2/2112/332 


= 2/2112/331 - 


2/1112/332 


= 2/3212/322 


— 2/2212/332 


= 2/3112/322 — yi\iy?,-ii = 


2/I11 


- 2/1112/331 


= 2/3112/321 - 


2/1112/332 


= 2/2222/321 


~ 2/2212/322 


= 2/2212/321 - 2/2112/322 = 


2/221 


— 2/2112/222 


= 2/2112/321 - 


2/1112/322 


= 2/2222/311 


- 2/2112/322 


= 2/2212/311 - 2/1112/322 = 


2/iii 


- 2/1112/221 


= 2/2112/311 - 


2/1112/321 


= 2/2112/221 


- 2/1112/222 


= 0. 



A simple verification yields that these equations define an irreducible variety of di- 
mension 3 and degree 9, which must be our eliminant. 



4.1.1 Coupling Constants 

In the above treatment, we have fixed a precise, known, form of the superpotential; 
essentially, we have fixed the coupling constants therein to be ea/37. Of the 27 possible 
gauge invariants which could contribute to W, we could take arbitrary complex con- 
stants to prefix each. Indeed, we could even have composite loops and thus infinite 
number of choices. Here, for convenience, we only consider minimal loops of length 3. 



We study the resulting moduli space of the elimination problem (4.2) by consider- 
ing different choices for the eap-y- First, we consider taking each e^^^ to be a random 
complex number. The numerical irreducible decomposition of the 9 generators involv- 
ing only the X variables yields six irreducible component of dimension 3, three sextic 



and three linear. Using (3.1), each component projects to a zero-dimensional variety, 
namely the origin. 

Next, we consider taking eaaa to be any nonzero complex number and all others to 
be zero. Then, the 9 generators involving only the X variables define a monomial ideal 
of dimension 3 that decomposes into 27 irreducible linear components. As above, each 
component projects to a zero-dimensional variety, namely the origin. 

Finally, if we take each e^fl^ to be 1, the 9 generators involving only the X variables 



decomposes into 3 linear components of dimension 7. Using (3.1), we know that each 
component projects to a 5 dimensional variety. Using the approach of § 3^, one 
readily obtains that each of these 5 dimensional varieties is generated by 15 linear and 
24 quadratic polynomials, and can verify that each has degree 12. The following are 
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the generators for one of these varieties with the generators for the other two obtained 
by a cychc permutation on the indices: 



yiii 


- 2/212 - 


2/213 


- 2/312 - 2/313 = 2/331 + 2/332 + 2/333 = 2/321 + 2/322 + 2/323 




yi2i 


— 2/222 — 


2/223 


- 2/322 - 2/323 = 2/311 ^ 


' 2/312 - 


^- 2/313 = 2/231 + 


2/232 - 


2/233 


— 


2/131 


- 2/232 ~ 


2/233 


- 2/332 - 2/333 = 2/221 ^ 


- 2/222 - 


^- 2/223 = 2/211 + 


2/212 - 


^ 2/213 






2/133 + 


2/233 


+ 2/333 = 2/132 ^ 


- 2/232 - 


^ 2/332 = 2/123 + 


2/223 - 


^ 2/323 






2/122 + 


2/222 


+ 2/322 = 2/113 ^ 


" 2/213 " 


^ 2/313 = 2/112 + 


2/212 - 


^2/312 


= 0, 


2/3232/332 — 


2/3222/333 




2/3132/332 — 2/3122/333 = 


2/2332/332 — 2/2322/333 = 


2/2232/332 — 


2/2222/333 = 


2/2132/332 — 


2/2122/333 




2/2332/323 — 2/2232/333 = 


2/2322/323 — 2/2222/333 = 


2/3132/322 — 


2/3122/323 = 


2/2332/322 — 


2/2222/333 




2/2322/322 — 2/2222/332 = 


2/2232/322 ^ 2/2222/323 = 


2/2132/322 — 


2/2122/323 = 


2/2332/313 — 


2/2132/333 




2/2322/313 — 2/2122/333 = 


2/2232/313 ~ 2/2132/323 = 


2/2222/313 — 


2/2122/323 = 


2/2332/312 — 


2/2122/333 




2/2322/312 — 2/2122/332 = 


2/2232/312 — 2/2122/323 = 


2/2222/312 — 


2/2122/322 = 


2/2132/312 - 


2/2122/313 




2/2232/232 — 2/2222/233 = 


2/2132/232 — 2/2122/233 = 


2/2132/222 — 


2/2122/223 = 



The Hilbert series for each irreducible component of the ehminant is 

1 + 7t + 4t2 

with the Hilbert series of their union being (1 + 15t + 18t^ + 2t^)/(l — t)^. 



4.2 Instantons on C 

The moduli space of instantons is an important but notoriously difficult object to study, 
for physicists and mathematicians alike. Even for flat (or, equivalently, the complex- 
ified C^) as the base space on which the instantons are realized as stable sheaves, the 



situation is unwieldy. Beautiful combinatorial approaches have been studied in 53 55 



Recently, some nice progress has been made in understanding the situation of /c = 1 



and 2 instantons with various gauge groups 56 , 57 , using Grobner basis techniques as 
well as direct integration using the Molien formula. In 33 , we alluded to some prelim- 
inary results on how one might approach the problem numerically. In this subsection, 
let us undertake a more systematic investigation. 



For concreteness, let us focus on A; G instantons with gauge group U (N). In 57 
it was realized that the instanton moduli space Mk,N is simply the Higgs branch of the 
vacuum moduli space of the following J\f = 1 quiver gauge theory in four dimensions. 
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whereby providing us the perfect formahsm as a testing ground: 



SU(N) 




N k 2 k 

i=l a, 6=1 a,/3=l a,b,c=l 



(4.3) 

In the above, we have fields which is a bi-fundamental of SU{N) x U{k), meaning 
that it is an X matrix of complex entries, likewise we have Qa? which is x A^. 
Moreover, we have {4>a=i,2)b which are adjoint fields under the U{k), meaing that they 
are k x k complex matrices and e"^ is the totally anti-symmetric matrix ^ ^ 
Finally, we also have $^ which is another adjoint of U{k). Throughout, we will use 
indices a,b,c, . . . = 1, . . . ,k to associate with U{k), i,j = 1, . . . , N to associate with 
SU{N), and a, /3 = 1,2. 

The moduli space in which we are interested is the Higgs branch wherein the vacuum 
expectation value of the field $ vanishes. Therefore, we have gauge invariant operators 
composed of the loops other than $ in the quiver diagram, which are of the following 
types: (1) Tr [(0i)^ ■ (02)'^] where ■ and Tr are k x k matrix multiplication and trace, 

N 

and the matrix exponents A, B = 0, . . . , k range so that < A + B < k; (2) ^ Q\Q\ 

i=l 

k 

for a, 6 = 1, . . . , /c; (3) E Q\{4>a)hQ] ior i, j = 1, . . . , N and a = 1, 2. Next, taking 

a,b=l 

the partial derivatives of W and setting $ to zero, the only non-trivial F-term is 

N 

E QaQ'i + 1^1^ ^2]a ioT a,b = 1, . . . ,k with the usual commutator forkxk matrices. 

i=l 

In summary, we have the following elimination problem. The indices are such that 
i, j = 1, . . . , N; a,b = 1, . . . , k; a = 1,2; and A, B = 0, . . . , k such that < A + B < k. 

Vab - Tr [(01 ) 

vl, - E Q\Q\ 
1 

k 



Given ideal 



v%. - E Q\{K)IQ) , 

a,h=\ 

N 
i=l 



1 \Va)b\ 



Eliminate : Q", 



(4.4) 
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Component-wise, we thus have 2k{N + k) variables coming from the quiver fields as 
well as 3k{k + l)/2 + 2N'^ coming from the y-variables. Hence we are performing 
elimination of 2k{N + k) variables in . 

The case {k,N) = (1,1) is an example that one can solve by hand. The one 
generator not involving the y variables is a monomial that defines two irreducible 
components, both of which project to the same plane in with Hilbert series 1 / (1— t)^. 

For {k,N) = (1,2), it is easy to verify that one generator not involving the y 
variables is irreducible. Thus, the eliminant is also irreducible that one can easily 
compute has dimension 4 and degree 6 with Hilbert series (1 + 4t + 1^)/ (1 — t)^. 

For {k,N) = (1,3), the one generator not involving the y variables is also irre- 
ducible. In this case, the eliminant is an irreducible variety of dimension 6 and degree 
30 with Hilbert series (1 + 12t + 15t^ + 2t^)/{l - t)^. 

For {k,N) = (2,1), the given ideal has three irreducible components, all of di- 
mension 8. However, the eliminant only has one irreducible component, which has 
dimension 6 and degree 6 with Hilbert series {1 + 2t + 2t^ + t^) /(I — t)^- 

For (A;, A^) = (2, 2), the four generators not involving the y variables define an irre- 



ducible variety of dimension 12 and degree 16. With (3.1), we have that the eliminant 



is irreducible of dimension 10. Using § 3.4, we find the eliminant has degree 72 with 
Hilbert series (1 + 4t + lOt^ + 19^3 + 2lt* + 12^^ + + t^)/{l - t)^^ . 

For {k,N) = (2,3), the four generators not involving the y variables define an 
irreducible variety of dimension 16 and degree 16 with the eliminant being irreducible 



of dimension 14. Using § 3.4, we find the eliminant has degree 2232 with Hilbert series 



(1 + lot + 55^2 + 183t3 + 410t^ + Qllt^ + 572t^ + 306t^ + 78t^ + Qt^)/{1 - ty\ 

For {k,N) = (3,1), the original ideal has four irreducible components, all of di- 
mension 15. However, the eliminant only has one irreducible component, which has 
dimension 10 and degree 180 with Hilbert series 



(1 + 7t + igt^ + 37t^ + 52t^ + 37t^ + 19t^ + It' + t^)/{l - t) 
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Finally, for (fc, A^) = (3, 2), the nine generators not involving the y variables define 
an irreducible variety of dimension 21 and degree 512. The eliminant is irreducible of 
dimension 18 and degree 1020 with Hilbert series 

(l+5t+15t^+34t^+62t^+96t^+130t^+157t^+170t*^+164t^+126t^°+52t^^+8t^^)/(l-t)^^ 
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4.3 Defining Algebraic Varieties 



It should be noted that our method of ehmination apphes to situations far beyond 
merely computing the moduli space of vacua for gauge theories; the setup is indeed 
applicable to problems in mathematics and physics alike. In this section, let us take a 
purely mathematical example which has over the years become important to the study 
of string compactifications. 

The earliest and still one of the most important dataset of Calabi-Yau three- 
folds, which constitute the corner-stone of compactifications of the ten-dimensional 
superstring, is the so-called complete intersection Calabi-Yau threefolds, or CICYs for 



short 58 , 59 . These manifolds are defined as smooth embeddings into products of 



complex projective spaces P"^ x . . . x P'^'= by polynomials of appropriate multi-degree; 

k 

the complete intersection criterion requires that there are ^ — 3 polynomials. We 

i=l 

will now take some explicit examples from this dataset and investigate some typical 
questions which arise. 



4.3.1 Embeddings into Projective Space 

Sometimes, we may wish to untangle the structure of the ambient space of product of 
projective spaces and embed our manifold into a single projective space. The advantage 
of doing so is that we have a single grading when dealing with bundles defined on the 
ambient projective space (which can then be restricted on the Calabi-Yau manifold) 
whereby facilitating computations of cohomology. The down-side is that generically 
we will lose the complete intersection property, though for purposes of computational 
algebraic geometry this is not particularly troublesome. 

Luckily, there is a standard technique of embedding product of projective spaces 
into a single one: the Segre embedding. To illustrate, let us consider the following 
famous Calabi-Yau manifold, the so-called bi-cubic (we take the form from [60, § 2.2]), 

2 2 

B{x,y) := xoXiX2^ yf + yoyiy2j2^i + ^oXiX2yoyiy2+ 

1=0 i=0 
2 2 

+ E ^hf+j + E E ^hi+kvf+jyi+j+k = 

i,j=0 i,j=Ok=±l 

(4.5) 

The matrix configuration specifies that we have a single polynomial of bi-degree (3, 3) 
in P^ X P^ whose respective projective coordinates we have labeled explicitly. We 



p2 o 

[xo:xi:X2] 
p2 o 



22 



have also given an example of this polynomial, as indicated. The subscripts on the 
coordinate variables are defined modulo 2. 

The Segre embedding here takes pair-wise products of the x, y projective coordi- 
nates, thus taking x into P*^: 

S : (xo,xi,a;2) x (2/0,1/1,2/2) — > Zm=o,...,& ■= (a;i?/j)ij=o,i,2 • (4.6) 

We wish to compute the defining equation B of the bi-cubic in P^. This, again, is 
described as an elimination problem: 



Given ideal 



zi - XoVi , 



- 3:22/2 
2 



> C C[a;o,i,2, 2/0,1,2; zo,. 



X0X1X2 yf + vovm + xoXiX2yoyiy2+ 



1=0 i=0 
2 2 



i,j=0 i,j=0 k=±l 



Eliminate 



3^0,1,2, 2/0,1,2 

Hence, we wish to eliminate 6 of the 15 variables. It is easy to verify that V{B) and 
the eliminant is irreducible of dimension 3. This implies that the eliminant is a cubic 
hypersurface inside the variety defined by the nine Pliicker relations: 

2:0^4 ~ ZiZ^ = ZqZ^ — Z2Z3 = ZqZ-j — ZiZq = ZqZs — Z2ZQ = ZiZ^ — Z2Z4 

= ZiZs — Z2Z7 = ZsZf — Z4ZQ = Z^Zs — Z^Zq = Z4ZS — Z5ZJ = 0. 

Since the Pliicker relations define a variety of dimension 4 and degree 6 in P^, the 



eliminant is a dimension 3 variety of degree 18 in P*^. The approach of § 3^ computes 
the final cubic relation: 

3 3 2 

Zq + ZqZiZ2 + ZqZiZ^ + 2:0^2^6 + ^0^3^6 + ^1 ~\~ ZiZq -\- Z\Z2Z4 + ZiZ^Z^ -\- Z1Z4ZJ 

+ ZiZq + ^2 + ^2^3 + ^2^7 + ^2^3 + ^2^4^5 + ^2^4^6 + ^2^5^8 + ^2^6^8 + ^2^7 

3 3 3 2 2 3 

+ ^3 + + Z^Z^Zq + + Z^Zk^Zj + Z^ZqZy + + Zc^Zq + Zk,Zq + Z^ZjZg + Zq 



+ ZqZjZq + Zj + Zg 



0. 



(4.7) 



4.3.2 Checking Smoothness 



A related problem to the above is whether one could efficiently ensure that our algebraic 
models for the manifolds are smooth. In the example of (4.5), we have chosen a 
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specific bi-cubic with fixed complex coefficients. In general, we can form ('^3^7^) ~ 
10 cubics in 3 variables; thus in general we could have 10 x 10 = 100 possible bi- 
cubic monomials, each of which, together prefixed by some arbitrary coefficient, could 
potentially contribute to out defining equation B{x,y). This is, of course, very much 



analogous to the situation considered in § 4.1.1 With an overall factor, we thus have 



a moduli space of projective dimension 99 possible manifolds: 

J2 Yl ^ijk,i'j'k' xix{x^y'oyiy^ = (4.8) 

< i, j, k < 3 < i' , j' , k' < 3 

1 + J + fc — 3 i' -\- j' -\- k' — 3 

We must impose conditions on the coefficients Aij^^i'j'k' such that the above represents 
a smooth threefold. Very often, we need to decide this efficiently and over large number 
of choices of the coefficients. Traditionally, this has been achieved in the literature by 
working over a coefficient field being ¥p for some prime p and checking over several 
small values of p. 

We will take a perhaps more gratifying approach and utilize our efficient numerical 



methods. Adhering to our example of (4.5) for concreteness, suppose we are given 



B{x, y; C, D) := xqXiX2 ^ yt + yoym + xoXiX2yoyiy2 + Cjxlyf^^ + 

j=0 i=0 ij=0 

2 

+ Dj^kxix^+kyl+^yi+j+k = (4.9) 

i,j=0 k=±l 

with 3 complex coefficients Co,i,2 and 6 complex coefficients -Do,i,2;±- Now we wish to 
check for what choices of these coefficients within our moduli space does the bi-cubic 
define a smooth Calabi-Yau threefold. Therefore, we need to check the gradient of 
partial derivatives with respect to the coordinates and see whether the simultaneous 
solution of these, together with B{x, y; C, D), indeed is the empty set in x P^. 

Once more, therefore, we are confronted with an elimination problem: 

r Bix,y;C,D) , | 
Given ideal : < ^^B{x,y;C,D) , > C C[xo,i,2, 2/0,1,2; Co,i,2, ^o,i,2;±] 
[ ^Bix,y;C,D) J 

Eliminate : Xo,i,2, 2/0,1,2 • (4.10) 



Following § 3^, we first compute a numerical irreducible decomposition for the 
given ideal. Over C, this yields 37 irreducible components of codimension 5, 36 of 
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which are hnear spaces and the other has degree 2101. Each component projects to a 
hypersurface in C^. For lu = (1 + iy^)/2 and a; = (1 — iy^)/2, each of the following 
twelve hyperplanes are the image of exactly 3 irreducible components: 

Co = 0, Ci = 0, C2 = 0, 

3(Co + C1 + C2 + Do+ + Do_ + Di+ + Di_ + D2+ + D^-) + 7 = 

3(Co + Ci + C2) - 3a;(L'o+ + ^1+ + ^2+) - 3a;(L>o- + + D2-) + 7 = 



3(Co + Ci + C2) - 3a;(L>o+ + ^1+ + ^2+) - 3c^(L>o- + ^1- + ^^2-) + 7 = 



3(Co + Ci + C2) + 3(L>o+ + i^o-) 


- 3a;(Di+ 


+ D2-) 


- 3a;(Di_ 


+ D2+) 
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= 



The other component projects to a degree 111 hypersurface in C^. Therefore, outside 
of these 13 hypersurfaces, the bi-cubic defines a smooth Calabi-Yau threefold. 

5 Conclusions and Outlook 

Computational algebraic geometry has been an immensely useful tool to study various 
theoretical physics problems, more prominently in string and gauge theories. However, 
the time has come when one needs to transition from the computation of relatively small 
models (say, number of variables and equations are of order 10) to the computation of 
more elaborate models. Unfortunately, the standard computational algebraic geometry 
methods run out of the steam in this transition because the Grobner basis techniques 
on which all these methods are based are exponential growth in running time and 
memory consumption. 

We have been introducing a novel numerical method called numerical algebraic ge- 
ometry (NAG), which efficiently resolves some shortcomings of the traditional Grobner 
basis technique including being completely parallelizable, to problems which involve 
large systems of polynomial equations, arising especially in theoretical physics. 

The typical question in gauge and string theory is concerned with the geometrical 
and particularly topological nature of the solution space, such as the dimension, the 
number of branches (irreducible components), how many isolated solutions, etc., for 
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this purpose, computational geometry is somewhat an over-kill since we usually assign 
random values to the parameters in any event and then perform Grobner reduction. 

For example, if we wish to check whether a geometrically engineered gauge theory 
indeed has a holographic dual being a certain Calabi-Yau manifold, we often set the 
coupling constants and mass parameters in the Lagrangian to generic values, compute 
the D- and F-flatness conditions, obtain a large set of polynomial equations, inter- 
pret it as an algebraic ideal, and then attempt to find topological quantities such as 
dimension and Hodge numbers. Therefore, if an efficient methodology were available 
to numerically resolve these issues then we can bypass the expensive Grobner basis 
calculation altogether. In this case, solutions can be computed to arbitrary accuracy 
which, in many situations, is sufficient for the application. However, if exact roots are 
indeed essential, one can always resort back to Grobner reductions or try exactness 
recovery methods such as [49]. 




The present paper is not only a part of our continuing program but it includes 
a new algorithm to compute a numerical irreducible decomposition of an eliminated 
ideal and (in the aCM case) also compute the important Hilbert series. We have also 
succinctly described how one of the most important problems in (supersymmetric) 
gauge theory, viz., finding the geometry of the vacuum moduli space, can be viewed as 
problems in elimination theory and applied our aforementioned algorithm to solve many 
examples ranging from geometrical engineering to instanton moduli spaces. We have 
even extended it to answering questions which is common-place is algebraic geometry: 
when is a given variety non-singular? 

The landscape of related open problems which confronts us is vast, and given our 
current tool we are optimistic in facing them. For example, the instanton moduli 
spaces for higher number of instanton number and for more complicated Lie groups 
are still mysterious. Calculating the decomposition and associated Hilbert series gives 
a handle on the partition function and counting. We have only demonstrated with some 
examples and marching toward these in general is an obvious direction. As another 
example, we have entered the age of string phenomenology where tremendously large 
databases of Calabi-Yau manifolds and related geometries are being created; issues 
such as checking smoothness and computing Hilbert functions are the very bread and 
butter of the field. Our numerical techniques should enormously help by bypassing the 
traditional Grobner reductions. As a final enticement, the geometry of the vacuum 
of the MSSM is yet unknown, and we will turn to this important computation in 
upcoming work. 
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